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We study preheating, i.e., non-perturbative resonant decay, of flat direction fields, con- 
- - - centratine; on MSSM flat directions and the right handed sneutrino. The diff'erence between 

cn; 

inflaton preheating and flaton preheating, is that the potential is more constraint in the 

o 

, latter case. The effects of a complex driving field, quartic couplings in the potential, and the 

presence of a thermal bath are important and cannot be neglected. 

Preheating of MSSM flat directions is typically delayed due to out-of-phase oscillations of 
^ the real and imaginary components and may be preceded by perturbative decay or Q-ball 
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X 



formation. Particle production due to the violation of adiabaticity is expected to be ineffi- 
cient due to back reaction effects. For a small initial sneutrino VEV, (N) < rriN/h with 



\l ' the mass of the right handed sneutrino and h a yakawa coupling, there are tachyonic insta- 

Tij" ' bilities. The Z?-term quartic couplings do not generate an effective mass for the tachyonic 

o : 

, modes, making it an efficient decay channel. It is unclear how thermal scattering affects the 

o 



resonance. 



PACS numbers: 98.80. Cq 



INTRODUCTION 



H ' During inflation the scale of quantum fluctuations is set by the Hubble constant. Scalar fields 

which are light with respect to the Hubble constant fluctuate freely, resulting in condensate forma- 
tion. Such condensates can play an important role in the evolution of the universe. A condensate 
can affect the thermal history of the universe, for example if its decay is accompanied by a large 
entropy production P|. The baryon asymmetry in the universe may originate from a conden- 
sate: Affleck-Dine baryogenesis utilizes a scalar condensate along a flat direction of the super- 
symmetric standard model P|, whereas Refs. discuss non-thermal leptogenesis from a decaying 
right-handed sneutrino condensate. Condensates can fragment into non-topological solitons, called 
Q-balls, which may have implications for dark matter Q, Oj. In the curvaton scenario, quantum 
fluctuations of a condensategenerate the density perturbations responsible for the observed CMB 
temperature anisotropy 0, UQ| . 

Field directions along which the effective mass vanishes, so-called flat directions, are a generic 
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feature of supersymmetric theories. There are many directions in field space along which the poten- 
tial vanishes identically at tree level. The classical degeneracy along the flat directions is protected 
from perturbative quantum corrections in the supersymmetric limit by the non-renormalization 
theorem P|. The flatness is lifted by soft terms from supersymmetry breaking. If the soft terms 
are sufficiently small during inflation the fields along the flat directions will be excited. Examples 
of flat directions in the MSSM are the HuL and udd directions ^ liol . Other examples of 
flat directions are the right-handed sneutrino and moduli fields. The masses of pseudo-Goldstone 
bosons are protected by an approximate global symmetry, and therefore can be kept naturally light 
during inflation. The Peccei-Quinn axion is an example of a pseudo-Goldstone boson that might 
condense during inflation. We will generically refer to flelds parameterizing flat directions as flaton 
fields, or simply flatons. 

Inflation erases all inhomogeneities along the flat direction, leaving only the zero-mode conden- 
sate. The vacuum expectation value (VEV) of the homogeneous mode can become large during 
inflation. In the post-inflationary epoch the field amplitude is initially damped by the expansion 
of the universe and remains essentially fixed. This stage ends when the Hubble constant becomes 
of the order of the fiaton mass, at which point the field starts oscillating around the minimum of 
the potential. This is similar to what happens in chaotic inflation models, where at the end of 
inflation the inflaton fleld starts oscillating around the minimum of the potential 11211 . 

Reheating of the universe through inflaton decay is a well studied problem fl3]. It is now 
understood that in many models non-perturbative processes, collectively known as preheating, can 
lead to rapid decay of the inflaton fleld 



141 . l^l . ^ , li , 3] • In preheating the decay of the inflaton 
occurs resonantly, leading to a rapid ampliflcation of one or more bosons to exponentially large 
occupation numbers. This process is eventually halted by the expansion of the universe, and by 
the back reaction of the produced quanta. The decay products scatter of the inflaton field leading 
to further decay of the oscillating inflaton zero-mode ^^,[2^. 

Flat direction condensates with large initial amplitudes may also undergo resonant decay. There 
is almost no mention of this possibility in the existing literature on cosmological scenarios based 
upon flat direction condensates. In this paper we analyze the possible occurrence of flaton pre- 
heating. We start our discussion in the next section with a review of preheating. We will highlight 
the differences between inflaton and flaton preheating. We restrict a detailed discussion of flaton 
preheating to two specific examples: Preheating of MSSM flat directions is the subject of section 
III, whereas section IV is devoted to resonant decay of the right-handed sneutrino. We conclude 
in section V. 
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II. PREHEATING 



Preheating with a real driving field has been discussed in detail in 
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Q]. Numerical stud- 
ies of preheating including the subsequent thermalization phase can be found in 

Q, Q- The 



generalization to complex fields is described in 



i3- 



We will summarize the main results. 



A. Real driving fields 

The flat direction condensate starts oscillating in its potential when the Hubble parameter 
becomes of the order of the flaton mass. Higher order terms in the potential rapidly become 
negligible, and the scalar potential can be written as 

with $ the field parameterizing the flat direction, and x some other scalar field it couples to. For 
the moment we take <^ real, postponing the discussion of preheating with a complex driving field 
to the next subsection. In supersymmetric theories the quartic coupling can arise from either a 
Yukawa coupling in the superpotential, or from gauge interactions. 

The equation of motion for the rescaled modes Xk = a in an expanding FRW universe 

is 

Xk+u;lXk = 0, (2) 

with k the comoving momentum, and a dot denotes differentiation with respect to time. Further, 
iol = k'^/a{tY + + h?^'^{t) + A, with A = -\{a/af' - |(a/a). The scale factor is a oc with 
p = 1/2 {p = 2/3) in a radiation (matter) dominated universe. Then A = |(^^ — + ^R- 

Soon after the onset of flaton oscillations A becomes negligible small. The flaton zero-mode is 
(J) = (j)cos{m^t) with the amplitude red shifting as oc a^^^^^. 



The mode equation can be brought in the form of a Mathieu equation 
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X'l + {Ak - 2q cos 2z)Xk = , (3) 



with z = m^t, and prime denotes differentiation with respect to z. Further 



fc^/a^ +rni 



X 

mi 



Am? 



(5) 
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An important feature of the solutions to the Mathieu equation is the existence of exponential in- 
stabilities. For q > 1, many resonance bands are excited. Preheating occurs in the broad resonance 
regime where particle production is efficient for modes with momenta k'^ < A — 2q. The occupation 
numbers of quantum fluctuations grows exponentially: oc exp(/ifcz) with significant exponent 
fik ~ 0.1. The typical momenta A:* of the particles produced is A;* = k^ax/V^ = n^-^Q^^'^ / V^, where 
we have assumed <^ ^m?^. 

In an expanding universe preheating attains a stochastic character. However the net result is 
still an exponential growth of ^-fluctuations. Preheating is halted by the expansion of the universe 
when q falls below unity. 

Effective preheating rests on two principles, violation of adiabaticity and Bose enhancement. 
Particle production results from non-adiabatic changes in the effective frequency of the Xk modes. 
Adiabaticity is violated when 

l^fcl ^ (6) 

which happens each time the flaton zero-mode goes through the minimum of the potential and 
changes rapidly. The occupation numbers of the decay quanta grow exponentially fast due to 
Bose-enhancement. As a result, preheating is robust. Resonant production occurs as long as the 
non-adiabaticity condition Eq. © is fulfilled. Adding additional fields or couplings (e.g. hm^x^) 
has little or no effect on the resonant period. However, the back reaction effects are very model 
dependent; when and how they become important depends on the specifics of the potential. 

B. Complex Driving fields 

Supersymmetric theories inherently involve complex scalar fields. Phase-dependent terms can 
arise naturally in the potential through soft SUSY-breaking terms. A relative phase between the 
oscillations of the real and imaginary components of the flaton fleld leads to a trajectory that is 
elliptic. The minimum amplitude is no longer |$|min 

= as it is for a real driving field, but instead 
I'l'lmin = b with b the semi- minor axis of the ellipse. This may prevent adiabaticity violation from 
occurring. 

We can decompose the complex driving field into real and imaginary components: <I> = <I>R-|-i<l>i. 
By a phase rotation the largest amplitude component of oscillation can be put in the real piece. 
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Then 

$1 = /0cos(m0t), (7) 

with / = b/a, the ratio of the semi-minor and semi-major axis of the elhptic trajectory. We will 
refer to / as the ellipticity of the orbit; note that / = \/l — with e the eccentricity of the 
ellipse. As before, the mode equation for the y-qu anta in the time-dependent <I>-background can 
be mapped into a Mathieu equation, with now 

AkU) = — - + 2qif), (8 

<p 

The (/-parameter for an elliptic trajectory is reduced by a factor (1 — f"^), with respect to the case 
of a pure oscillatory trajectory. 

Adiabaticity violation as defined in Eq. (jB)) occurs when 

^+ml + h'f4>' + (1 - f)h''^l < ((1 - f)h'<^nM^f' , (10) 

The ellipticity is negligible small when for typical momenta k ^ k^t, the h? f'^cj)^ term in the above 
equation can be neglected. This is the case for 

f^fR= 2g(o)i/4 ' 

and preheating proceeds as for a real driving field. For larger ellipticities the term proportional 
to in Eq. dominates. This leads to an upper bound on q{0) for which broad resonance is 
effective: 

qiO) <qc^ . (12) 



The reason is that for larger values of g(0) the semi- minor axis of the ellipse is large, b oc f\/q{0), 
and adiabaticity violation does not occur. For large initial q- values resonant decay is delayed until 
the expansion of the universe red shifts the q-parameter below the critical value. As discussed 
below, preheating can be very efficient if at its onset q(f) > 10'^, which requires / < 0.1. Note that 
for large ellipticities / = 1 — e > 0.5 the upper bound requires q ^ e, and preheating never takes 
place. 
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In supersymmetric theories bosonic preheating is generically accompanied by fermionic preheat- 
ing, since a Yukawa couphng in the superpotential leads also to a fermion coupling of the form 



£ 9 (mf + h(l))iljil), (13) 

with mf the fermion mass. Resonant production of fermions has been studied in As long as 
the flaton amplitude is larger than ruf/h (i.e., g ^ 1), and the ellipticity is sufficiently small, there 
is an instant during each flaton oscillation that the effective mass of the fermion vanishes, and 
fermions are produced. Within about ten oscillations the fermion occupation number is saturated 
at a time-averaged value ~ 1/2 for momenta within a Fermi sphere of radius fcf ~ m^q^^'^. The 
Pauli-exclusion principle forbids occupation numbers beyond one. The back reaction of fermions 



can catalyze bosonic preheating 



1^. 



Another interesting property of supersymmetric theories with regard to preheating was observed 



m 



12|, who studied a superpotential of the form W = ^S"^ + ^im^ + h<f>)x^- After diagonalizing 
the mass matrix for the real and imaginary component of Xi it turns out that one of the masses 
can become tachyonic during part of the oscillation period of (p. As a result, quantum fluctuations 
of Xk grow exponentially. We will refer to this kind of instability as tachyonic preheating. The 
occurrence of tachyonic modes is quite model dependent. 



C. Back reaction effects 



Which back reaction effects play an important role, and at what stage, is rather model depen- 
dent. We will list here four different back reaction effects for bosonic preheating. 

• When the quantum fluctuations of x grow exponentially large, the effective flaton mass 
"^eff ~ "''0 + 5'^(x^) becomes dominated by the variance term. The oscillation frequency of $ 
rapidly increases, and energy is rapidly dumped into x-pai'ticles until (7 ~ 1 and preheating is 
halted by the expansion of the universe. At the end of this phase (x^) ~ 0^, and occupation 
numbers are enormously large ~ 10^ h~^. Initial values q > 10^ are needed for this rapid 
energy transfer to occur; for smaller values preheating is halted by the expansion of the 
universe before the variance term comes to dominate the flaton mass [isl ]. 

• The non-zero variance (x^) can also induce an effective mass for the x-field itself if quartic 
couplings Ax'^ are present in the potential. Preheating of particles with masses m^^ ^ 



7 



is generically inefficient. Tfiis may fialt prefieating before the stage of rapid energy transfer, 
rendering preheating inefficient. 

• The decay quanta scatter with the zero mode, leading to exponential amphfication of flaton 
quanta. However, scattering only becomes important when the occupation numbers of the Xk 
quanta become exponentially large, in the very last stages of efficient preheating. Scattering 
leads to further decay of the zero-mode. 

• Flaton preheating can occur in the presence of a thermal bath. If the interaction rate for 
scattering of % quanta with the particles in the thermal bath is much greater than the flaton 
oscillation frequency, resonance modes are depopulated rapidly on the relevant time scale. 
There is no Bose enhancement, and resonant production is inefficient. 

D. Flaton vs. Inflaton 

In this subsection we consider the differences between resonant decay of generic flat direction 
fields and the usual considered case of resonant decay of the inflaton field. 

• The energy density stored in the inflaton field dominates the energy density in the universe 
during the epoch of preheating. As a result, the expansion rate of the universe is set by 
the inflaton field itself, and is initially that of matter domination. The flaton fleld on the 
other hand is generically subdominant at decay. It evolves in a fixed background, which can 
be either radiation or matter dominated. This difference in universe evolution leads to only 
small changes in the various quantities, leaving the order of magnitude estimates unchanged. 

• The universe is reheated by the decay of the inflaton field. Thus inflaton preheating occurs in 
an empty universe. This is not necessarily the case for flaton preheating. Note in this respect 
that even before the reheating process of inflation is completed there is a dilute plasma with 
temperature T ~ (T^MpiH)^/'^, where Tr is the reheat temperature of the universe. The 
effects of the thermal bath should be taken into account for flaton preheating. 

• The inflaton is usually chosen to be a gauge singlet of the standard model, its properties 
largely unconstrained by particle physics or experimental data. Moreover, preheating is 
mostly studied in the simplest setting: the inflaton is a real field, with only a quartic inter- 
action term. Flat direction fields naturally appear in SUSY theories. This means that for 
flaton preheating we can no longer ignore the complexity of the driving field. Furthermore, 
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we restrict our attention to two specific examples, namely the right handed sneutrino and 
MSSM flat directions. For these flatons the potential is constrained by supersymmetry and 
particle physics, and there is no freedom in ignoring/ adding unwanted/wanted terms. 



III. MSSM FLAT DIRECTIONS 



Directions in field space along which the potential vanishes identically are called flat directions. 
In supersymmetric theories there are generically many flat directions at the classical level. The 
classical degeneracy along flat directions is protected from perturbative quantum corrections in the 
supersymmetric limit by the non-renormalization theorem. The flatness is lifted by soft terms from 
supersymmetry breaking and non-renormalizable operators. During inflation the inflaton potential 
dominates the energy density in the universe. The non-zero energy density breaks supersymmetry 
leading to soft masses and ^-terms for scalar fields 22] ■ In this section we will discuss MSSM flat 



directions in detail, whereas the right handed sneutrino field is the subject of the next section. 

The MSSM flat directions can be parametrized by gauge invariant operators, X, formed from 
the product of m chiral superfields making up the flat direction. Defining X = cc/)™", the effective 
potential for the MSSM flat direction is of the form 



V = + cnH^m^ + ;^„_3 Aa>» + h.c.J + ' '^,1, . (14) 

Here M is the cutoff scale, typically the GUT or the Planck scale. In the MSSM with parity 
conservation, most flat direction are lifted by n = 4, 5 or 6 non-renormalizable operators. The 
flattest one is lifted by n = 9. The Hubble induced terms proportional to H are the soft terms from 
SUSY breaking by the finite energy density in the universe. Vm and the term proportional to the 
gravitino mass are the MSSM soft mass and A-term respectively. The mass term depends on 
the SUSY breaking scheme. For gravity mediated SUSY breaking Vm = my^\^\'^ at tree level. In 
the case of gauge mediation Vm = m^log(l + |<I>p/m^), where ~ 1 — 100 TeV. 

Condensation along the flat direction occurs during inflation, if (1) the flaton mass is much 
smaller than the Hubble constant and quantum fluctuations can lead to large initial field values; 
or if (2) the effective flaton mass squared is negative and the field settles at a minimum away 
from the origin. The natural size for the soft masses during inflation is of the order of the Hubble 
constant P,13|- Non-minimal Kahler potentials can induce negative soft (mass)^ terms, with 
ch < 0, realizing possibility (2). Scenario (1) can be realized in the context of e.g. D-term 
inflation [2^ or no-scale type supergravity models where symmetries forbid soft mass terms 



FIG. 1: The cUipticity / as a function of Oa — da for different Hubble induced masses m? = cH^ . Solid lines 
correspond to c = —5, dashed lines to c = —1, short dashed lines to c = —.25 and dotted lines to c = 10^^. 
All plots is for n = 4 in a radiation dominated background. 

at tree level and c//, a <C 1 naturally. 

Writing <I> = -^(pe^^, the 6 dependent part of the potential reads 

Va-^ ^^{\a\Hcosi9a + ex + ne) + \A\m3/2Cos{eA + ex + ne)) , (15) 

with 6a, Oa and 6x the phase of a, A and A respectively. During inflation the 0-field will settle 
in one of the minima of the o-term if a ^ 0(1)- For a <^ 1 the radial motion is dominated by 
quantum fluctuations. 

In the post-inflationary epoch the evolution of the flaton field is given by its equations of motion: 
(j) + J>H(l) + dV/d4> = 0. The (j) field is damped as long as H > negative (mass)^ {ch < 0) 

the damping is critical, and the field closely tracks its instantaneous minimum 

'^min(i^)-( ^™"'' )'^" (16) 

with P some numerical constant depending on a, ch and n. Quantum fluctuations of fields with 
a positive (mass)^ {ch > 0) saturate at V{cl)) ~ H'^ 25[, and initial amplitudes (j)o ~ Hj/m 



are expected, with Hj the Hubble constant during inflation. The field is over damped at the 
end of inflation, and remains essentially fixed until 0o 

^ 0min; and it starts slow rolling m the 

non-perturbative potential. 

Eventually the Hubble induced mass becomes equal to the soft mass H ~ \/V^i, and the flaton 
starts to oscillate in the potential well. This is the moment ellipticity is created. In gravity 
mediated SUSY breaking this phase starts when H ~ 7713/2- At this time the Hubble induced 
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FIG. 2: EUipticity for F-terni inflation a 7^ and n = 4 (solid lines); Z?-terni inflation a = and ri = 4 
(dashed lines); radiation domination and n = 6 (short dashed lines); and for matter domination and n = 6 
(dotted lines). For all plots c = 0. 

and MSSM ^-terms are of comparable magnitude and there is a torque in the angular direction if 
7^ ^a- We expect ellipticities of order unity. In the gauge mediated case 7713/2 < GeV and the 
MSSM A-term is small compared to the Hubble induced one at the onset of oscillations. Therefore 
s.^le. .lM.cities ..e expected. The eOiptieit—tote, to a constant when H ^ 0.1™, 

This is confirmed by numerical calculations [2a| . In gravity mediated SUSY breaking the ellip- 
ticity is / < 0.5 — 0.2 for n = 4 — 6. In more than half of the parameter space {9a — 9a) G [0, 27r] 
the ellipticity is / < 0.1. The parameters chosen in these simulations are M = Mpi and 
1^1 = \a\ = —\ch\ = 1- The results are independent of the gravitino mass and the cutoff scale. 
The ellipticities obtained in the gauge mediated case are smaller / < 0.1 — 0.01 for n = 4 — 6. 
These results are for 771,3/2 = 10~^, Id^'^m^; for much larger values the behavior as in the gravity 
mediated case. 

We have extended the numerical calculations of [2^ to study the elliptical trajectory for more 
general parameters, focusing on a mass term of the form Vm = ^^I^P- We found that the ellipticity 
has the same order of magnitude independent of the sign of the Hubble induced masses, i.e., it is 
independent of ch, as shown in Fig. 1. Further, the same results are obtained in both F-term and 
D-term inflation, see Fig. 2. In the latter no Hubble induced terms are generated, and a = 0. 

Ellipticities do depend sensitively on the ratio of the gravitino to flaton mass, and also on the 
initial amplitude of the flaton field. The ellipticity is suppressed by a factor 7^3/2/7770. This is as 
expected, as in the limit 7773/2/777,^ — > no potential in the angular direction is generated. This 
is shown in Fig. 3. For negative (mass)^ the field is trapped in the minimum of the potential, 
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FIG. 3: EUipticity for different ratios of the gravitino to flaton mass. Solid lines correspond to m^/^/m — 3, 
dashed lines to m^jilni = 1, and dotted lines to = 0.1. For all plots c = 0. 

given by Eq. (|l(ij) . But for positive (mass)'^ it is possible to have an initial amplitude differing from 
(pmin- Ellipticities decrease rapidly for (pQ < (pmim as shown in Fig. 4. For the parameter values 



7713/2 ~ TeV and M = Mpi, (p^^i^i 
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10^^ GeV for n = 4 — 6 at the onset of oscillations. 



Quantum fluctuations saturate at a field value lower than 0min for Hubble constants during inflation 
Hr < 10^ - 10^ GeV forn = 4 - 6. 

The flaton has D-teim couplings to the fields making up the polynomial X, which are orhtogonal 



to the direction which aquires the VEV |2|^. In addition, the flaton has F-term couplings to fields 
not present in X. The interaction term is of the form h?(f)'^x^, with h either a gauge or yukawa 
coupling. The q-parameter is q{f) = (1 — /^)/i^i?!>^/(4m^). The gauge and Yukawa couplings in the 
MSSM range between /i ~ 1 — 10^®, m^^ = 7713/2 ~ ITeV, and the field amplitude at if ~ is 
given by Eq. (fTC)) . Typically g(0) » 1. We can then distinguish the following possibilities. 

1. The ellipticity is negligible small. Preheating starts when H ~ and can be effective, 
depending on the back reaction effects. Such small ellipticities are obtained for fine-tuned 
parameters Oa ^ da, for (po <^ (/'min as is possible for positive (mass)^ and low scale inflation, 
or in the absence of ^-terms. Generically, supergravity corrections will induce non-zero 
^-terms. 

2. The ellipticity is small / < 0.1. This is the case in gauge mediated SUSY breaking, and also 
in most of the parameter space for gravity mediation. Moreover, ellipticities are suppressed 
if < 1 and/or 0o/0min < 1- Preheating starts off in the broad-band regime. 
However, its onset is delayed until q{0) ^ Qc ^ 1/(16/^), and therefore might not occur if 
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the condensate decays perturbatively, via thermal scattering or through fragmentation into 
Q-balls before this time. Note that the decay quanta that will be first exited are the fields 
corresponding to the smallest g-parameter, i.e., the fields with the smallest coupling to the 
flaton. 

Denote the decay width of the condensate by ~ /^m^. For perturbative decay [5 = h? 
with h a Yukawa or gauge coupling. For temperatures higher than the effective mass of the 
particles coupling to the flaton, i.e., T > [m-^cS ~ ^4*^ there is a thermal bath of x-particles. 
Thermal scattering can lead to decay of the condensate; in this case (3 ~ h?a{T /m^) with 
a = (7^/47r. Note that in the presence of a thermal bath the flaton mass should be corrected 
by the thermal contribution in all formulas Sm^ ~ h?T'^. If it dominates, the thermal mass 
induces early oscillations. Now the field amplitude, and therefore not only decreases 
because of the red shift, but also because the effective mass decreases. This will speed up 
the onset of preheating, since g(0) < qc earlier. Thermal effects are particulary important 
for n = 4 directions with small Yukawa couplings h < 10~^ 27]. Since the thermal history is 
rather model dependent, we will not pursue this issue further. Finally, if the flaton potential 
grows less than (jp' the condensate is unstable against fragmentation into Q-balls. Numerical 
simulations indicate that Q-ball formation takes place at i7 ~ 10~^ — lQ~^m^ in gravity and 
gauge mediated scenarios respectively P|. 

If we denote the Hubble constant at the moment (7(0) = qc with He, then preheating can play 
a role if ^ F^. In a matter dominated universe this translates into a^^^f'^h < 'm^/(j)mm, 
with (prain evaluated at -ff ~ m^. The universe is matter (radiation) dominated before 
(after) inflaton decay. For example, for ~ TeV and M ~ Mpi preheating occurs before 
perturbative decay if h^^f^ < 10^^, (10^^^) for n = 4, (6). Preheating precedes Q-ball 
formation if f^h < 10~^ — 10"^ for n = 4 — 6 and gravity mediated SUSY breaking, and 
f'^h < 10~^ — 10^^ for n = 4 — 6 and gauge mediated SUSY breaking. 

3. The ellipticity is appreciable 0.2 < / < 1. Preheating is delayed. If it takes place at all, 
it is shut off by the expansion of the universe before the back reaction becomes important. 
Large, order one, ellipticities are possible in gravity mediated SUSY breaking with n = 4 
non-renormalizable operators. 
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FIG. 4: EUipticity for different values of the initial amplitude. Solid lines correspond to (/)o/0min = 1, dashed 
lines to ^o/'/'min — 5, short dashed lines to ^o/'/'min = 0.1, and dotted lines to (f>Q/(l)min = 0.05. For all plots 
c = 0. 

A. backreaction 

Quartic self-couplings of the decay mode V 3 h?x^ will induce an effective mass term my. ~ 
h?{'X^) as the variance grows. All MSSM bosons (except for the right handed sneutrino) are charged 
under the standard model gauge group, and have a quartic interaction with h the gauge coupling. In 
addition, the superpotential term W 3 hcpx^ leads to a quartic interaction in the .F-term potential. 
It has been found that for a potential of the form V{(l)) + m^x^ + h?(fp"x^ resonant production is 
only effective for couplings h? > 10~'^(mv/m^)^ [l^ . Therefore, bosonic preheating shuts off when 
("^x)efT ~ "^(/>- From = g4>o{x'^) [l^l, it follows ~ m^(j)o/h ^ p,^ ~ '"^^'/'o- Thus at the of 
preheating, only a small fraction of the energy stored in the flaton zero mode will be transferred 
to the decay products. 

No effective mass term for the fermion superpartners of x generated. However, the Pauli 
exclusion principle forbids occupation numbers exceeding unity. The typical energy transferred to 
the fermions is kp ~ /i^m^0g, which is only significant for large couplings. 

The typical interaction rate for MSSM particles in a thermal bath is of the order a ~ a^T. 
The relevant time scale in preheating is . It follows that thermal scattering is important for 

> m^/y^HMpi for T < Tr, and > m^/{T^HMpiy/^ for T > Tr, with Tr the reheat 
temperature of the universe. Effective scattering kills Bose-Enhancement, but can increase the 
efficiency of fermionic preheating, since in between each moment of particle production the decay 
products are scattered out of the resonance bands. 
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IV. RIGHT HANDED SNEUTRINO FIELD 

The right-handed neutrino and its scalar partner appear in grand unified theories (GUT) based 
upon SO{W). At first sight it appears the sneutrino direction is hfted by a D-term potential due 
to the 50(10) gauge interactions, preventing a large VEV for the sneutrino field A^. The crucial 
observation made in _3'] however is that the D-term decouples from the potential of N as long as 
the value of N is smaller than the 50(10) breaking scale. The sneutrino field condenses during 
inflation if <^ H^, with a maximum amplitude set by the GUT scale. 

The superpotential containing the three right-handed sneutrino fields (after B — L breaking) is 

3 3 
^ = Y. l^MiNiNi + ^ Ko^NiL^H , (17) 

i=l i,a=l 

here Aj are the three sneutrino superfields, La the three left-handed lepton doublets, H is shorthand 
for Hu the down Higgs superfield, and hia are the Yukawa couplings. We start our discussion 
of preheating by considering a single generation. Writing the scalar potential in terms of the 
component fields = L") and H'^ = {H+ 

Vf = m%N^ + \hf\N\WLf + \H\^) + h*mNN{L°H° - H+f + c.c. (18) 

= ^(|l0|2 _ |^0|2)2 ^ 9^^l0.j^~ ^ H+*H^f + c.c, (19) 
8 2 

Vs = mL\L\'' + mH\H\^ + (Bml/^ + hH^)N^ + ^l!^-t^A" + c.c. (20) 

with |Lp = |L°p + and \H\'^ = \H'^\^ + Vp and Vd are the F and L>-term potential 

respectively, and Vs is the soft SUSY breaking potential. The coefficients B,b, A,a ~ 0(1), 
gw = e/ sm6\Y and gz = ^/ sin 9w cos Ow- The a, 6-terms are from SUSY breaking by the finite 
energy density in the universe, whereas the A, i3-terms arise from low energy SUSY breaking. A 
is the ultraviolet cutoff scale, typically A ~ Mqut or A ~ Mpi. If R-parity is conserved the lowest 
order ^-term for the sneutrino has n = 4. Furthermore, we have assumed that non-renormalizable 
operators are sub-dominant. An inclusion of these operators will not affect the results in any 
essential way. 

Since the i?-terms do not contain a coupling between the real and imaginary parts of A, it 
will not induce an ellipticity. The sole source of ellipticity are the A-terms. During infiation the 
Hubble induced ^-term is dominant and the the 9 field quickly settles in one of the minima of 
cos{9a + 9\ + n9). At the onset of oscillations, when H mj\f, the low energy A-term is negligible 
small for gravitino masses 771,3/2 <^ m]\r, and the field remains essentially fixed in its 0- minimum. 
The field amplitude red shift with time and all ^-terms quickly become sub-dominant; hereafter 
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no further ellipticity is generated. We expect therefore that for sneutrino masses rriN ^ "^3/2 the 
eUipticity is neghgible small, and we can use the theory of preheating for a real driving field. This 
is confirmed by our numerical calculations which give / ^ 0. 1 (771,3/2 /jn- at), see Fig. 4. Therefore we 
take real. For now we also take the yukawa coupling h real, deferring a discussion of a complex 
coupling to section IIVBI 

We decompose the slepton and Higgs fields in its real and imaginary component. Diagonalizing 



the mass matrix gives the following mass eigenvalues and eigenstates: 

M^^ =hN{hN ±mN), fovXR±=L%±H% (21) 

Mf^ = hN{hN ±mN), for X/± = L? =F , (22) 

= hN{hN ± niN), for Xr± = L^^ H+, (23) 

= hN{hN ± mTv), for Xi± = LJ ± H+ , (24) 



where the unbarred (barred) quantities correspond to the neutral (charged) fields. 

Due to both the linear and quadratic term in N we cannot map the mode equation for the 
decay products in terms of a Mathieu equation. However, adiabaticity violation and thus particle 
production occurs when goes through zero, and 

\^k\>ujl (25) 

During most of the resonance time the quadratic term dominates over the linear term. Hence, 
the physics is well captured by this term alone, and the problem can be rewritten in terms of a 
Mathieu equation with q = h? Nq / {Amj^) . Particle production occurs in a time interval 5tnon-ad ~ 

The (mass)^ eigenvalues are not positive definite. During part of the A^-oscillation the masses 
can go negative, and quantum fluctuations of the X_ fields grow exponentially due to the 
tachyonic instability. For KNq < rrin or equivalently (7^1, the various mass terms become 
tachyonic during half of the sneutrino oscillation, and tachyonic preheating is very efficient. On 
the other hand, for 9 ^ 1 the mass terms only become tachyonic during the small time interval 
^haxAi ~ 1/(^-^0)- Since ((5t)non-ad ^ ((^Otach) it is expected that in this regime tachyonic particle 
production is sub-dominant. Particles with masses up to m^/4 can be produced during tachyonic 
preheating. . 

The (7-parameter for sneutrino decay is = h'^ Nq / Amj^ . In the seesaw mechanism the neutrino 
masses are related to the Yukawa couplings and sneutrino masses through rriu = ^ , where 
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{Hy} = 174 GeV x sin/3 and tan/3 = {Hu)/{Hd}. The solar and atmospheric neutrinos have masses 
in the range 10""^ — 10~^ eV. Typical values of the q-parameter are 

Typically g 3> 1, and wc expect non-adiabatic preheating to happen generically. Note however 
that there is no lower bound on the lightest neutrino mass, and therefore some of the Yukawa 
couplings might be arbitrarily small. Moreover it might be that N gives a negligible contribution 
to the neutrino mass, which are dominated by their couplings to the the sneutrinos in the other 
two families. In those cases the g-parameter can be lowered to arbitrary low values. 

A. Back reaction effects 

The quartic interaction terms can be written in terms of the mass eigenstates: 

£,2 ^2 _ _ _ _ 

Vd = -^{Xi+Xj- + Xr+Xr-) — ^{Xr+Xi- + Xii-Xi+ — Xi+Xr- — Xi-Xr+) 

+^{Xi+Xi_ + Xi_Xj+ + Xr+Xr_ + XR_XR+f (27) 

/,2 

(28) 

Consider first the case g S> 1. Then tachyonic instabilities are ineffective, and the eigenstates Xi 
all grow with the same rate. The D-term potential generates effective mass terms for all modes 
{m'x)cS ~ g'^{X'^). Just as in the MSSM case, this effective mass will halt preheating before 
a significant amount of energy is transferred from the sneutrino zero mode to the sleptons and 
Higgses. This conclusion can be avoided if one of the X-fields has a large VEV before the start 
of preheating. For example, if {Xj^} is initially large, it gives a large mass to the and X/_ 
fields. As a result, these eigenstates are not exited during preheating, and no effective mass for the 
field is generated. 

The case < 1 is quite different. Tachyonic preheating is effective, and only the tachyonic 
eigenstates X_ are produced. The D-term will not generate an effective mass state for these fields. 

The _F-tcrm potential does lead to an effective mass term mx_ h'^{X'^). Preheating halts when 
the effective mass term becomes of the order of the right handed sneutrino mass. Before this time 
a significant amount of energy can be transferred from the sneutrino zero mode to the tachyonic 
fields. 
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The last back reaction effect to be considered is thermal scattering. The typical interaction rate 
for Higgs, selectron and left handed sneutrino particles in a thermal bath is of the order T ~ T, F ~ 
a^T and T ~ G\T^ respectively, with Gp the fermi constant. The relevant time scale in preheating 
is m^^. It follows that scattering is important, both for the Higgs and slepton fields, even before 
inflaton decay when there is a dilute plasma with temperature T ~ (T^ffMpi)^/^. Decay quanta 
are scattered rapidly out of the resonance bands. It requires further studies to determine how 
effective thermal scattering is in killing the resononance. The reason is that tachyonic preheating 
can be very rapidly itself: In the absence of a thermal bath the zero mode typically decays within 
one or two oscillations I2SII. 



B. CP violation 



To study CP violation during preheating we allow the coupling constant h to be complex. 
Particle production due to violation of non-adiabaticity is dominated by the quartic CP conserving 
term h?N'^'x^. Hence we expect CP violation to be small. 

To study CP violation for the tachyonic decay mode, we decompose h = hji+ihi and diagonalize 
the mass matrix. The mass eigenvalues and eigenstates are: 

Ml± = \h\N{\h\N ± M) for Xr^ = L%± j^/?" ^ (29) 

Mf^ = \h\Ni\h\N ± M) for = L? ^ j^F? ^ ^LhI (30) 

(31) 

Similar expression holds for the charged lepton and charged higgs fields. It follows that tachyonic 
preheating will not generate a lepton asymmetry. 

Finally we consider the fermionic decay modes. The mass matrix for the neutrino and higgsino, 
in the basis {v^ H^)-, is of the form: 

/ n hNX 

(32) 

Although the fermion mass matrix is hermitian, CP is broken since for a complex coupling since 
M 7^ M* . A lepton asymmetry is expected in fermionic preheating, as discussed in reference 
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V. DISCUSSION 

We studied non-perturbative, resonant decay of flat direction fields, concentrating on MSSM flat 
directions and the right handed sneutrino. The difference between inflaton preheating and flaton 
preheating, is that the potential is more constraint in the latter case. The effects of a complex 
driving field, and of quartic couplings in the potential are important and cannot be neglected. 
Moreover, preheating occurs in the presence of a thermal bath. Effective scattering can kill Bose- 
enhancement effects. 

Ellipticities of MSSM flat directions condensate are generically of order one, thereby delaying 
the onset of preheating. Preheating may be preceded by perturbative or thermal decay, or by 
the formation of Q-balls. If preheating does occur it is generically ineffective due to the quartic 
self-couplings of the decay modes, and due to the presence of the thermal bath. 

The ellipticity for the right handed sneutrino field is negligible small. Particle production due 
to the violation of adiabaticity is generically expected to occur. However, just as in the MSSM 
case, it is expected to be inefficient due to back reaction effects. The new feature in the sneutrino 
potential is that for small values q = h^N^ / {Amiq) < 1, there are tachyonic instabilities. The 
L'-term quartic couplings do not generate an effective mass for the tachyonic modes, making it an 
efficient decay channel. It is unclear how thermal scattering affects the resonance. 
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